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INEQUALITIES BETWEEN SIZE AND CHARGE FOR BODIES AND THE 
EXISTENCE OF BLACK HOLES DUE TO CONCENTRATION OF CHARGE 

MARCUS A. KHURI 


Abstract. A universal inequality that bounds the charge of a body by its size is presented, and is 
proven as a consequence of the Einstein equations in the context of initial data sets which satisfy 
an appropriate energy condition. We also present a general sufficient condition for the formation of 
black holes due to concentration of charge, and discuss the physical relevance of these results. 


1. Introduction 

It is well known that black holes of a fixed size can only support a certain amount of charge, 
depending on the horizon aresQ [8l ll4l[T8] . Here we propose a similar statement for arbitrary charged 
bodies which do not lie inside a black hole. Namely, let H be a compact spacelike hypersurface in a 
spacetime which satisfies a suitable energy condition. If H lies in the domain of outer communication 
then there exists a universal constant C such that 

(1.1) Charge(H) < C ■ Size(n), 

where a precise definition of size will be given later. Thus all bodies, from elementary particles to 
astronomical objects, can only support a fixed amount of charge depending on their size, or rather 
they must be sufficiently large depending on their charge. Similar results and inequalities have 
recently been obtained [Zmsille] where the role of charge is replaced by angular momentum, that is 

(1.2) AM(0) < C ■ Size(H) 

if H is not inside a black hole. The constant C in (|l.ll) will be a multiple of c^jy/G, where c is the 
speed of light and G is the gravitational constant, and Size(fl) will be measured in units of length. 

In |19j . it was shown that if the amount of angular momentum of a body sufficiently exceeds 
its size, then the body must be contained in a black hole. In this paper, we will also establish 
such a criterion for black hole existence focusing instead on the role of charge. More precisely, 
if the opposite inequality of o holds, then Q must be contained in a black hole. It follows 
that concentration of charge alone can result in gravitational collapse. This statement is naturally 
motivated by intuition, since large amounts of charge are associated with strong electromagnetic 
fields, and high concentration of matter fields is known to result in black hole formation. This last 
statement is referred to as the Hoop Conjecture [3T|, and is related to the Trapped Surface Conjecture 
m- These conjectures have received much attention, although they are not fully resolved. Many 
works [DiaEllISKITl [2011211133] involve strong hypotheses including spherical symmetry or maximal 
slices, and some results [ini|29l[35] are not useful for initial data with low extrinsic curvature. 


The author acknowledges the support of NSF Grant DMS-1308753. 

^This refers to apparent horizons. The same statement is not necessarily true for event horizons |22| . 
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2. Explicit Formulation 


Let (M, g, k, E, B) be initial data for the Einstein-Maxwell equations which includes a Riemannian 
3-manifold M with complete metric g, symmetric 2-tensor k denoting extrinsic curvature, and vector 
fields E and B representing the electric and magnetic fields. It is assumed that the stress-energy 
tensor is decomposed into two parts, one for the electromagnetic field, Tem, 

and the other for the remaining matter fields, Tm- If is the timelike unit normal to the slice, 
then /i = T’^^na'fib and J* = T^^'Ua represent the energy and momentum densities of all matter 
fields, whereas g,M = T^Uarib and J\j = T^Ua give the energy and momentum densities of the 
non-electromagnetic matter fields. The initial data must satisfy the constraints 


(2.1) 

Jm = J -\- -—E X R, 
dvr 

where 


(2.2) 

^4 k- R+{Trgkf l/cp 

J = div{k {Trgk)g), 


with R the scalar curvature of g, and {E x B)i = eijiE^B^ the cross product; here e is the volume 
form of g. Recall also that the electric and magnetic fields are obtained from the field strength by 
Ei = Fin and Bi = —^eijiF^K The following inequality will be referred to as the charged dominant 
energy condition 


(2.3) 


k-M > \Jm\- 


Let n be a body, that is, a connected open subset of M with compact closure and smooth boundary 
dkl. The sum of the squares of its electric and magnetic charges yields the square of total charge, 
which (using Gaussian units) is given by 


(2.4) 




divEdkjjf 


— I divBdujg 


in 


We now describe measurements of the size of the body 14. In this regard, two definitions of 
radius will be important. Namely, in [29] Schoen and Yau defined a radius TIsy{^), that may be 
characterized as the largest (minor) radius among all tori that can be embedded in fl, and which 
was central to their condition for the existence of black holes due to compression of matter. In more 
detail, if T is a simple closed curve that encloses a disk in H, and r is the greatest distance from 
T with the property that all points within this distance combine to form an embedded torus in fl, 
then TZsy{^) is equal to the maximum r from any such curve T. Another important radius, the O 
Murchadha radius TloM{kl), is defined [24| as the radius of the largest stable minimal surface that 
can be embedded in H. Here, radius of the surface means the largest distance from a point in the 
surface to the boundary dQ, as measured by the induced metric on the surface. This is formulated 
most simply when dQ is mean convex (has positive mean curvature), so that geometric measure 
theory guarantees the existence of many smooth least area surfaces contained in H. In general, the 
O Murchadha radius gives a larger measure of size than the Schoen/Yau radius 


( 2 . 5 ) 


Eom{^) > Esy{^)- 
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Both radii measure well the size of a ball of radius a in flat space. Namely, for this body TZsy = o/2 
and TZqm = o- However, their measurement for a torus of major radius a and minor radius b is less 
accurate: TZsy = b/2, TZqm = b. In particular, for a torus, this measurement of size does not take 
into account the major radius. This leads to a problem if one tries to establish an inequality of the 
form with the notion of size given in terms of either of these radii. For instance, in the weak 

field limit a torus of large major radius a but small minor radius b could still support a large amount 
of charge, since its surface area and volume may be large, while the measure of its size in terms of 
the radii is small. For this reason, we choose a notion of size which incorporates surface area |9H| as 
well. That is, in the precise version of inequality dni, size is defined by 

( 2 . 6 ) SiMQ) = 

where TZ{Q) represents either the Schoen/Yau radius or the 6 Murchadha radius. Lastly, it should 
be mentioned that the radius TZqm gives an accurate measurement for highly dense spherical bodies 
[3 [21]. Thus even in a strong gravitational field, this measurement is on the order of the area radius. 
Nevertheless it would be desirable to eventually replace the role of the Schoen/Yau and O Murchadha 
radii with a more simple quantity. However our results here depend on theorems proved for these 
more complicated measurements, and it is not clear how to proceed without them. 

With these notions of charge and size, we may formulate a rigorous description of inequality (jl.ip . 
except for the constant C which will be given in the next section. In order to articulate the black 
hole existence result, we must replace event horizons with the quasilocal notion of apparent horizons. 
This is due to the fact that event horizons cannot be located in initial data without knowledge of the 
full spacetime development, whereas apparent horizons may be identified directly from the initial 
data. Let S' C M be a 2-surface and let := Hs ± Trgk be the null expansions, which assess 
potency of the gravitational field; here Hs denotes mean curvature with regard to the unit outward 
normal. Geometrically, the null expansions arise from the first variation of area in the outward future 
(0+) and outward past (0_) null directions, and as such they are the rate of change of area for a 
shell of radiation discharged by S in these directions. If 0+ < 0 {6- < 0), then S is called a future 
(past) trapped surface, and the gravitational field is considered to be strong in this vicinity. The 
boundaries of future (past) trapped regions are referred to as future (past) apparent horizons, and 
solve the equation 0+ = 0 (0_ = 0). Cosmic censorship implies that apparent horizons are typically 
contained inside black holes [32], and thus in many situations they may be used in place of event 
horizons. We will show that if the opposite inequality of (|l.ip is valid, then an apparent horizon 
must exist within the initial data. 


3. Inequalities Relating Size and Charge of Bodies 

In this section, inequalities of the form (jl.ip will be established, both in the maximal case {Trgk = 
0) and in the general case. The inequality obtained in the maximal case is stronger, as the universal 
coefficient C in this case is smaller. However the inequality obtained for general initial data will be 
used to obtain the criterion for black hole existence, described in the next section. We begin with an 
important observation which will be used in both cases. Let H be a body as described in the previous 
section, then the total charge may be estimated in terms of the energy and momentum densities as 
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follows 


Q^ = 


1 

47r 


E ■ vdOa I + 


'90 


1 

47r 


B ■ vda. 


'90 


(3.1) 


=1^^ / [(l-®P + “ Stt//+ SttI Jm|) + 87r(/i - I JmD] do-g 


990 

where in the last step the charged dominant energy condition (12.3p was used, and u is the unit outer 
normal to <90. 


Theorem 3.1. Let (M, k, E, B) he a maximal {Tvgk = 0) initial data set for the Einstein-Maxwell 
equations. Then for any body Ll C M with constant energy density g and satisfying the charged 
dominant energy condition (I23D, the following inequality holds 


(3.2) 


IQI < 


\dn\ 


12G7^(^1)’ 


where TZ{Ll) denotes the Schoen/Yau radius TZsy{^)- Moreover, if in addition the boundary dLl is 
mean convex, then TZ{Ll) denotes the O Murchadha radius TZom{M). 


Proof. It suffices to estimate the integral on the right-hand side of (j3.ip . 

(3.3) [ ig-\JM\)dag< [ gdag = g\dLl\. 

Jan Jan 

In light of the maximal assumption and the constancy of g, Theorem 1 
yield 


We have 


of [29] may be applied to 


(3.4) 

It follows from (13.111 that 


^ - 6G7^5Y(^1)2■ 


(3.5) 


< 




Now consider the case when the boundary dLl is mean convex. It was pointed out in m, that 
under this additional hypothesis, the estimate (j3.4p holds with the O Murchadha radius 

4 

TTC 

It follows that (|3.2I1 holds with the O Murchadha radius. Note that (|2.5p implies that this is a better 
result than the estimate with the Schoen/Yau radius. □ 


We will now obtain an inequality between the size and charge of bodies without the maximal 
assumption. Here we will employ a technique developed by Schoen and Yau in [281129] . which reduces 
certain problems for general initial data back to the maximal setting. The idea is that in the maximal 
setting, nonnegative scalar curvature i? > 0 is guaranteed from the dominant energy condition, and 
it is this nonnegativity which is fundamental for establishing many geometric inequalities, such as 
the positive mass theorem or (13.411 in the proof of Theorem 13.11 Thus, it is natural to deform the 
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initial data metric <7 to a new unphysical metric g whose scalar curvature satisfies i? > 0 , at least 
in a weak sense. This is accomplished in [2H] by setting = gij + VifVjf■, which is the induced 
metric on the graph t = f{x) in the 4-dimensional product manifold M x M, where / satisfies the so 
called Jang equation 


(3.7) 



rp \ ( 

1 + lv/IV IvVi + lv/P 



= 0 , 


with P = g^^Vjf. The purpose of this equation is to guarantee that the scalar curvature of g is 
weakly nonnegative, in fact it is given by the following formula O [ 6 l [28] 

(3.8) R = - J (v)) + \h-k\‘^+ 2\q\l - 2diVg{q). 

Here divg is the divergence operator, h is the second fundamental form of the graph, and 


(3.9) 


fi 

Vi + lv/P’ 




If the dominant energy condition g > \J\ is valid, then each term on the right-hand side of (|3.8I) 
is clearly nonnegative, except perhaps the divergence term; hence we may view R as being weakly 
nonnegative, which is sufficient for most purposes. 

Let us assume that the Jang equation (I3.7|l possesses a smooth solution in H. Measurement of 
the accumulation of matter fields, which is needed to estimate the right-hand side of (|3.ip . may be 
obtained by assessing the concentration of scalar curvature for the unphysical metric 'g. This in turn 
may be accomplished by estimating the principal Dirichlet eigenvalue 


(3.10) 


Ai — 


In 


of the operator — ^R, where (j) is the principal eigenfunction. From the weak nonnegativity of 
the scalar curvature, we may integrate by parts and use the two nonnegative terms |V(/>p and 
to find 


(3.11) 


^ SttG , . ^ 

Ai > -^nnn(/r - |J|) =: A. 


Here we have also used the fact that |u| < 1, so that g — J{v) > g — \J\. 

The estimate (|3.1ip for the principal eigenvalue allows an application of Proposition 1 of |29] (if 
A 74 0), from which we obtain 


(3.12) 


Rsvi^) < 



where the radius TZsy pertains to g. Since the metric g measures lengths larger than does the metric 
g, it follows that TZsy > TZsy- Furthermore, let V’ £ C°°(H) be an arbitrary positive function, then 
dividing and multiplying A“^ by f^ipdujg [f^ig — iJDipdujg) produces 

A-i < c^Cq 4 H^g 

~ SttG f^(g - |J|)V’dwg’ 


(3.13) 
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where Cq = when ^ — \ J\ > 0 in O, and Cq = oo when // — | J| = 0 somewhere in Q. 


Therefore 

(3.14) 


(// - I JDV’dWg < 


Sttc^Co !^^dug 


16G TZsvi^y' 

and this leads to a general correlation between the size and charge of bodies. Before stating this 
result, we record integral forms of the charged dominant energy condition and an ‘enhanced’ energy 
condition on the boundary 


(3.15) 


/ ifJ-M -\JM\)dag >0, / inM -\J\)dag >0, 

JdQ Jdn 


and define a constant Cf = Co where J is Jm or J depending on whether the first or 

second condition of (I3.15jl is satisfied, respectively. 

Theorem 3.2. Let (M, g,k, E, B) be an initial data set for the Einstein-Maxwell equations, which 
contains no compact apparent horizons. Assume that either M is asymptotically flat, or has a strongly 
untrapped boundary, that is Hqm > \TrQM^\- Then for any body Q C M satisfying the dominant 
energy condition /r > | J|, and one of the two energy conditions ()3.15l) on dLl, the following inequality 
holds 

‘ \dn\ 


(3.16) 


IQI <Ci 


3c4 


32G 7^5y(^l) 

Proof. The assumptions concerning the asymptotics of M or its boundary imply that a strongly 
untrapped 2-surface is present in the initial data, and hence the Dirichlet problem [29] for the Jang 
equation admits a solution with / = 0 on such a surface. In addition, / must be a smooth solution 
due to the lack of apparent horizons. The arguments above now apply, so that (I3.14p is valid. 

Let us now assume that the first energy condition of (j3.15D holds. Then (|3.ip is valid and we find 
that 

\dO.\ 


Q" < 


27r 


Ian 


ih - \JM\)dcrg 


^ ^ lanih - \ JM\)dag ^ maxan(/i - \ Jm\) |5n| 


(3.17) 

Furthermore by choosing 

(3.18) 
we obtain 

19 ^''^ 77sy(0)2 

from ()3.14p . Together, (I3.17P and (13.191) yield the stated conclusion. 

Alternatively, if the second energy condition of (13.151) holds, then add and subtract | J| instead of 
\Jm\ in (13.ip to obtain 

\dn\ 


fnih \T\)d^g 


/ { h -\ JM \) d ( Tg < 

Jan 


mino(^ - I J|) |fl| 

\dn\ 


Sttc^Ci 


(3.20) 

By choosing 

(3.21) 


Q" < 


27r 


[ {h-\J\)d(Tg. 
Jan 


if = 


_ fanih - \J\)d(yg ^ maxaQ(/i - | J|) |5n| 


in^h-\J\)du:g minn(^-|J|) |SJ| 


a similar argument yields (13.161) . 


□ 
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Remark 3.3. The first energy condition in (I3.15|) is satisfactory, as it is a weaker version of the well 
known energy condition (j2.3p (used for instance in the positive mass theorem with charge) which 
states that the non-electromagnetic matter fields satisfy the dominant energy condition. On the 
other hand, the second energy condition in (I3.15P is not well motivated, and we view it as a technical 
assumption that could perhaps be removed with further investigation. It should also be pointed out 
that the constant in (I3.16P is independent of the particular matter model, as long as it satisfies the 
appropriate energy condition. 

This result generalizes Theorem 13.11 which involves strong hypotheses such as the assumption of 
maximal data and constant matter density. Note, however, that the inequality of Theorem 13.21 is 

weaker than that of Theorem 13.11 since the constant Ci appearing in (j3.16p is generally larger 

than the constant of (j3.2p . Thus, although two undesirable hypotheses have been removed, 

the resulting inequality is not optimal; the problem of finding the optimal constant is currently being 
investigated by Dain et al. [9] in the context of spherical symmetry. We also mention that the 
charged bodies constructed by Bonnor in [3j, all of which satisfy (|3.16l) . could potentially be useful 
in this pursuit. It turns out that the difference in the constants just described is related to a black 
hole existence result which we now explain. 


4. Criteria for Black Hole Formation 

The dependence of Theorem 13.21 on solutions of the Jang equation ()3.7p inherently produces a 
black hole existence result. This is due to the fact that solutions are regular except possibly at 
apparent horizons, where the graph t = f{x) tends to blow-up in the form of a cylinder (see |15] . 
[28j). In other words, if it can be shown that the Jang equation does not possess a regular solution, 
then an apparent horizon must be present in the initial data. This method for producing black holes 
was initially used by Schoen and Yau in [29]. Here we will use it to obtain a criterion for black hole 
existence due to concentration of charge. 


Theorem 4.1. Let {M, g,k, E, B) be an initial data set for the Einstein-Maxwell equations, such 
that either M is asymptotically flat, or has a strongly untrapped boundary, that is Hqm > \TrQMk\. 
//H C M is a body satisfying the dominant energy condition g > \J\ and one of the two energy 
conditions (|3.15p . with 


(4.1) 


IQI >Ci 


3c4 \dn\ 


32G TZsYi^)' 

then M contains an apparent horizon of spherical topology which encloses a region that intersects H. 


Proof. As in the proof of Theorem 13.21 the assumptions on the boundary of M or its asymptotics 
imply the existence of a solution to the Dirichlet boundary value problem for the Jang equation, 
with / = 0 on dM or on an appropriate coordinate sphere in the asymptotic end. If the solution 
were regular, then by Theorem 13.21 the opposite inequality of ()4.ip would hold. Since this is not the 
case, we conclude that the solution is not regular, and hence the existence of an apparent horizon is 
guaranteed. Moreover, among apparent horizons arising from the blow-up of Jang’s equation, there 
is at least one (which is outermost) with spherical topology |28j . If the region which is enclosed by 
this apparent horizon does not intersect H, then the solution of Jang’s equation is smooth over H. 
This would then imply by the proof of Theorem 13.21 that (I3.16P holds, yielding a contradiction. □ 
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Whether concentrated charge leading to gravitational collapse is a naturally occurring phenome¬ 
non, seems to be an intriguing open question. In the next section we will comment on some physical 
aspects of this problem. As for the theoretical part, it should be pointed out that the proposed 
criteria is not satisfied in the maximal case, since as mentioned previously, the universal constant 
in inequality (13.2p is smaller than the constant in (13.161) and (|4.ip . A similar relation between the 
maximal and nonmaximal cases holds with regards to Schoen and Yau’s condition for black hole cre¬ 
ation. Therefore, sufficient amounts of extrinsic curvature are needed for the Schoen/Yau condition, 
as well as the hypotheses of Theorem 14.11 to be fulfilled. 

With this intuition, we now show how to construct examples of initial data satisfying the conditions 
for black hole existence. Fix an asymptotically flat metric gf on M ~ M^, and set 11 = i?i(0) 
to be the unit ball. Let ( 61 , 62 , 63 ) be an orthonormal frame and set kij := k{ei,ej) = 0 for all 
ihj) 7 ^ (1)1)) (2) 2), and ku = k 22 = I3‘^ for some parameter /3 > 0. It follows that 

(4.2) ^ = ^(ij + 2A:nfc22)~/3") | J| ~/3'. 

For large /3 this yields Ci ~ 1, and also hm > |^| if we choose E,B ^ /3. It may be arranged so 
that divE = divB = 47 r /3 in 11, which implies that \Q\ = ay/2\kl\. Since the right-hand side of 
(j4.ip is independent of /3, it follows that ()4.ip is satisfied for large j3. Lastly, by extending k, E, 
and B outside of H to be asymptotically flat, all the hypotheses of Theorem 14.11 are satisfied and 
an apparent horizon must be present. Since the constant 2/3^ represents Tvgk, we conclude that 
large traces of extrinsic curvature facilitate the formation of apparent horizons in this setting. Note 
that there is quite a bit of freedom in this construction, since for instance the metric g is essentially 
arbitrary. Of course, this is an abstract and ad hoc procedure which is not of physical interest, but 
it shows that there are many configurations satisfying the criteria of this result. It should also be 
pointed out that besides the Maxwell field, the matter models present in this construction are not 
given explicitly. However, it seems that a combination of dust (or perfect fluid) with a charged scalar 
field might fit. Finally we mention that since g — \J\ ~ /3^, and the radius TZsy{^) is fixed, this 
data also satisfies the hypotheses of the Schoen/Yau black hole existence criterion [29]. Moreover a 
similar construction in axisymmetry, with the role of charge replaced by angular momentum, yields 
examples of data satisfying the criteria of m- 


5. Physical Relevance 

The inequalities between size and charge for bodies, as well as the black hole existence criterion 
proven above, are predictions of Einstein’s theory and hence should be contrasted with observational 
evidence and other theories. Let us consider bodies which are approximately spherical in shape, so 
that the ratio of boundary area to radius is on the order of the radius TZ. Then in general terms, 
what we have shown is that for stable bodies 

and that if the opposite inequality holds then the body should undergo gravitational collapse. Here 
< should be interpreted in terms of order of magnitude, and ke ^ 9 x 10^Nm'^C~‘^ is Coulomb’s 
constant so that (15.ip is expressed in SI units, as opposed to Gaussian units used in previous sections. 
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Consider now an electron. It has a classical radius of TZe ~ 2.8 x 10 Moreover, since 

G « 6.67 X Nm'^kg~'^ and c « 3 x it follows that 

c2 

(5.2) ~ lOOC. 

\ rv^ijr 

Therefore, since the charge of an electron \Qf,\ ~ 1.6 x lO'^^C, we hnd that (|5.ip is satisfied. 

According to the principle of charge quantization, the charge of a body is an integer multiple of 
the elementary charge (charge of an electron). Thus, \Qe\ is the smallest amount of charge that a 
body can possess. Using this fact in (15.Ih . we find that the classical theory imposes the following 
minimum size for a body 

(5.3) 77-0 = -^—^\Qe\ ~ 1-4 X 10“^®m, 

which is on the order of the Planck length Ip = (^) ~ 1.6 x It then appears to be a 

remarkable self consistency of the Einstein field equations that they predict a minimum length on 
the order of magnitude of the Planck length, if we assume the principle of charge quantization. 

On the other hand, we may consider bodies of astronomical scale such as stars. The study of the 
effects of electric charge in isolated gravitating systems goes back to Rosseland [27] and Eddington 
m]. It was shown that since electrons are rather less massive than protons, electrons tend to escape 
more frequently, as part of the solar wind. This induces a net positive charge in the star, which then 
yields an attractive force on electrons trying to escape. Eventually an equilibrium of these forces is 
established, resulting in a net positive charge on the order of ~ 1OO(M/M0)C' [16] . where M is the 
mass of the star. Thus, for typical stars, net charge is sufficiently small to be considered insignificant, 
and they certainly satisfy inequality ()5.ip . However, as pointed out by Witten [3l|, it is theoretically 
possible to have stars made of absolutely stable strange quark matter. These are highly dense bodies, 
which have masses and radii similar to those of neutron stars. They are also capable of possessing 
large amounts of charge |23|, and thus are candidates to violate ()5.ip . Consider such a star with 
charge \Q\ = lO^'^C' and radius TZ = 10^ m as considered in [23|. We have 

so that dsn) is still satisfied, although it is nearly violated. Moreover, the black hole existence 
criterion associated with (15.ip asserts that a star of this radius can only support a charge of \Q\ ~ 
beyond which the system will collapse to form a black hole; this is consistent with the findings 
of [25], obtained numerically with different methods. Lastly, we mention that magnetic charge is 
also included on the left hand side of and thus it would be interesting to contrast the above 

results with empirical evidence associated with magnetic charge. 
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